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ALFVEN WAVES IN SPIRAL INTERPLANETARY FIELD 

t 

Y.C. Whang 
Abstract 


This paper presents a theoretical study of the Alfven 
waves in the spiral interplanetary magnetic field. The 
Alfven waves under consideration are arbitrary, large 
amplitude, non-monochromatic , microscale waves of any 
polarization. They superpose on a mesoscale background 
flow of thermally anisotropic plasma. Using WKB approxi- 
mation, an analytical solution for the amplitude vectors 
is obtained as a function of the background flow properties 
density, velocity, Alfven speed, thermal anisotropy, and 
the spiral angle. The necessary condition for the validity 
of the WKB solution is discussed. The intensity of fluctua- 
tions is calculated as a function of heliocentric distance. 
The relative intensity of fluctuations compared with the 
magnitude of the background field has Its maximum in the 
region near 1 a.u. Thus outside of this region, the solar 
wind is less turbulent. Due to attenuation of microscale 
Alfven waves, fluctuation energy is converted into the kine 
tic energy of the solar wind. Its contribution amounts to 
about 20 per cent of the kinetic energy at 1 a.u. 


I. 


INTRODUCTION 


Power spectra and discontinuities of the interplanetary 
magnetic field have been studied by Burlaga (1968, 1969), 
Coleman (1966),, Sari and Ness (1969), Siscoe et al. (1968), 
and others. These studies reveal that a considerable amount 
Of microscale fluctuations is present in the interplanetary 
medium. 

The frequency dependence of the power spectrum based on 
interplanetary magnetic field data were performed' to 
study the statistical properties of the magnetic field 
fluctuations. Their results show that the mean amplitude 
of the fluctuations is usually a substantial fraction of 
the background field magnitude. The power densities of 
field fluctuations at a given frequency for the disturbed 
periods vary from 2 to 15 times those in the quiet periods. 
Except on very quiet days, the field fluctuations are not 
isotropic. The power levels in the field magnitude gene- 
rally are substantially lower than those of the components, 
and the fluctuations are more intense in the direction 
transverse to the average background field. 

Hydromagnetic discontinuities in the interplanetary 
medium are frequently observed. Most discontinuities 
are directional discontinuities ( rotational or tangential 
discontinuities). Studying the spectral slope of the 
Pioneer 6 magnetic field data. Sari and Ness (1969) found that 
at times the power spectrum in the microscale range is 
dominated by the presence of directional discontinuities 
in the microstructure of the interplanetary magnetic field. 
One of the directional discontinuities, the rotational 
discontinuities are also known as Alfven shocks or trans- 
verse shocks (Colburn and Sonett, 1966). Like smooth 
Alfven waves, the rotational discontinuities propagate 
along the field with Alfven speed. Fluctuations in the 
velocity and field vectors for rotational discontinuities 
and for smoother Alfven waves are also governed by the 
same dynamical relations. Thus in a broader sense, Alfven 
waves include rotational discontinuities. 


2 


Recently Belcher et al. (1969), and Belcher and Davis ( 1971 ) 
studied the interplanetary fluctuations in velocity and 
magnetic field vectors using the simultaneous plasma and 
field data from Marinen 5. They showed that Alfven waves 
with periods in the range of 

T 4 

1 ''“'10-10 seconds 

propagate outward from the sun, and they are the major 
source of microscale fluctuations in the interplanetary 
medium. 

This paper presents a theoretical study on the propagation 
of Alfven waves in the interplanetary space. The Alfven 
waves under study are of arbitrary form, large-amplitude 
and non-monochromatic waves, include smooth waves as 
well as rotational discontinuities , and are not restricted 
to any special polarizations. We consider that the micro- 
scale Alfven waves with periods in the range of 10 to 10 
seconds superpose on a smooth background interplanetary 
field and plasma flow. In the time scale of —-10 seconds, 
rotational discontinuities actually appear like continuous 
current sheets. The background interplanetary magnetic 
field has a spiral configuration with its radial component 
inversely proportional to the square of heliocentric distance 
The thermal anisotropy of the solar-wind plasma is included 
in the theory. V/e consider that the temporal and the 
spatial variations of the background field and plasma 
flow are measured by a mesoscale rime of 

t >, 10^ seconds 

2 

or its equivalent length scale of 
L >. 4xl0 12 cm 

near 1 a.u. Thus we are dealing with a theoretical problem 

involving two time scales. The important concept of the 

various interplanetary time scales was first introduced by 

Burlaga and Ness (1968). (They divided the microscale and 

3 4 

the mesoscale at ''*4x10 seconds instead of ~10 seconds.) 
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We obtain analytical solution for the mesoscale attenuation 
of the Alfven waves in the interplanetary space. It includes 
the effect of the spiral angle <*>. Some quantitative dis- 
cussion of the result is also given. 


II WKB APPROXIMATION 


Let V and B respectively denote the momentary value of 
the' velocity and the magnetic field vector, and U and E^, 
the mesoscale average value of V and B with respect to 
time at a point fixed in space 


U = < V > and B = <B > 


Here the angle brackets are defined as 


<X(r,tX>: 


t{_ 


■T/2 


T/2 


X(r,t+t) dr 


For practical reason, we take the averaging time T to be 
a finite time interval. This interval must be at least 
several times the upper limit of the microscale period of 
the wave motion. Now we may write 

V = U + u Cl) 

and 


B = B 0 + b 


( 2 ) 


where u and b are microscale fluctuation vectors. The 


governing equations for the variation of u. and b are 
' 3u 


at + ^ B o .?b = 

3b 

— ,+ U .7b - B Q . Vu = -(V- U)b 




(3) 

(4) 


Two kinds of terms are contained in the above two 
equations, the time-derivative of the fluctuation vectors 
and the derivatives of various quantities in the direction 
of the mesoscale average vectors U and B . The divergence 
of U represents the fractional change in mass density fol- 
lowing the mesoscale motion of the flow 

yU = - -(|- + u.y)p 

p 3 1 — 


V 


The effect of plasma thermal anisotropy is included in 
Equation (3) through the presence of a parameter 5 which 
is defined, as 

5 = 1- (Pjj-Pjl) /B 2 = 1 - 0(T|j-T_ L )/2T (5) 

The parameters C and 0 which characterize the thermal 

properties of the magnetized plasma are calculated based 
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on the time-average value of B , i.e. 

<B 2 > = B 2 +<b 2 > # ( 6 ) 

When the two terms on the right hand side of Equation (3) 
and (4) are both zero, the equations describe the pro- 
pagation of hydromagnetic waves along the convecting 
magnetic field lines with Alfven speed 

A = (5/4irp) ' B q (7) 

and the two fluctuation vectors are related by 

u = HF(5/4 ir p) 1/2 b. (8) 

In this case. Equations (3) and (4) are reduced to 

^ +. (U±Aej).v) (b,u) = 0 (9) 

where ej = j^/B^ In ( 8 ) and (9) the upper sign refers 
to anti-parallel disturbances and the lower one to parallel 
disturbances. When an arbitrary disturbance is convected 
and propagated as Alfven waves, b and u remain as two 
constant vectors. In other words, the amplitude of Alfven 
waves remains unchanged in a region where 5 and p 
are uniform. 

The characteristic length for the terms on the left hand 
side of (3) and (4) is the wave length of a monochromatic 
Alfven wave, X. Let L be the scale height for the 
variation of the background flow. Then the relative im- 
portance of the terms on the right hand side of (3) and 
(4) compared with those on the left hand side may be 
measured by a parameter 

y = X/L 
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y is a small parameter (y<l) when we consider that the 

of 

wave length is^microscale and the scale height L is of 
mesoscale. The small parameter y is equivalent to 
when the background, flqw, is also uns.teady in mesoscale. 

When o<y<l, the microscale Alfven waves are modified 
due to the slow variation of the background flow in the 
mesoscale. We may write each fluctuation vector for the 
solution of (3) and (4) as the product of a mesoscale 
amplitude vector and a microscale wave function f. Using 
WKB approximation the amplitude vectors are expanded as 
series in power of y with the leading terms of the series 
satisfying the zeroth order solution (8). The microscale 
wave function f represents an arbitrary function of its 
argument satisfying the equation for traveling Alfven 
waves in a general way as though the mesoscale effect 
were absent, 

).vjf = 0 (10) 

The spatial and the temporal variation of f are charac- 
terized by microscales, X and A special case for 
monochromatic sinusoidal waves, 

f = exp(iS - iujt) 

with 

- w - (_U ± Ae^)» vS 

Under these approximations, we may write the two fluctuation 
/vectors as 


b = (bj+yb2+»»*«)f 

(11) 

and 


u = (n 1 +pU 2 + • • . . )f 

(12) 

where 


Uq= ? (^/Uttp 

(13) 


and f remains as an arbitrary function. 

The spatial and the temporal variation of b^'s and u^'s are 
characterized by mesoscales L and T 2 * 


-ft + (H 4 A ^i 




Substituting (10)- (13) into (3), and Keeping terms to 
first order in p , we obtain 


6 


pA^ J. H2±— 2y— 1’ 7 .^' + ‘ f ("|^ + .y *'7-±Ae^ * v )b^_ + 

+ fbijj^(|t + M* v± 2 A £i‘ 7) ^ - i^ ( |t + y #v) pj= ° (i4) 

Similarly, we obtain from (4) 

-yA^-J— ^ u 2 ± -2)-l* 7f + f( ft + ^* 7 1 a ®i ,7 ^i + 

+ if ± 7^~r' 7 S “ “(ft + y* 7 i f ®i* 7 ^pj = 0 (15) 

Add (14) to (15) and rearrange the terms, we obtain a 
singel differential equation for the variation of b-j_ 


.1/2 1/2,9 


1/4 -3/4, 


( |_ + y.7)( 5 i,H p ° /H b x )± Ae x .v(c 


3/4 -1/4 


b x ) = 


= 0 


( 16 ) 


*ti 


This equation holds in any coordinate system and with no 
restriction on the models of the background flow. No 
assumption has been made concerning the polarization of 
the fluctuation vectors. 

In the region where the background flow has a very high 
Alfven number (U>>A), Equation (l£>) gives 


1/4 -3/4 

K p b 1 « constant 


( 17 ) 


On the other limit, when the flow is steady and the Alfven 
speed is very large (A>>U), the variation of wave amplitudes 
is governed by 


3/4 -1/4, 
£ P 




constant 


( 18 ) 
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III PROPAGATION IN SPIRAL INTERPLANETARY FIELD 

The solar-wind flow has a very lov/ Alfven number in the 

corona. The Alfven number of the solar wind increases 

* 

to approximate 10 near the orbit of the earth, it continues 
to increase beyond 1 a.u. We would like to find an exact 
solution of Equation (16) to cover the whole range of 
.Alfven numbers for the variation of the amplitude vectors 
of the Alfven waves superposed on the solar wind. 

Consider a steady, spherically-symmetrical model of the 
solar wind, the background field has a spiral configuration 
In general, the radial flow velocity U is not parallel to 
the direction of the spiral field. Let $ denote the angle 
between U and e^. Then we can write (16) in the following 
convenient form 


_ A ,d ,. 3/4 -1/4, . 1/4 -1/4. .d ,.1/2 1/2., 

(U±A cos4>)^(£ p - b^) - £ p U^;( 5 p 


= 0 


(19) 


From the equation of divergence B, the radial component 
of the background magnetic field, B Q cos*, is proportional 
to r"" 2 , we can write 

j^(p/c) 1/2 Ar 2 cos* 




( 20 ) 


The equation of continuity is 

|j(pUr 2 ) = 0 

From (20) and (21) we obtain 


( 21 ) 


.-1/2 -1/2 d ,.1/2 1/2,, A cos* d , U 
? p a ?' 5 p >+ D - 1 d? ( A cos ~» ) = 0 (22) 

ubstituting (22) i nto (19), we obtain 


, U . ,,d ,.3/4 -1/4^ , 
( A “ESi? * n d? <5 » ^i> 


3/4 -1/ 4d . U 

^ dr a cos^- 1~ 


0 


which can be integrated to give 


t 3/4p- 1 /4 ( _Jj__ ± 1)£i = constant (23) 

Here the amplitude vector is given as a function of £,p, 

U, A and the spiral angle $. Since no limitation on the 
component of the background field transverse to the radial 
direction has been made in the above derivation, the result 
holds for very general spiral fields. 

Two special cases of Equation (23) have been studied in 
literatures. Farraro and Plumpton (1961) derived a solution 
for the vertical propagation of Alfven waves in a horizontal- 
ly stratified, stationary atmosphere. They assumed that 
the fluctuation vectors are linearly polarized in a hori- 
zontal direction and the medium is thermally isotropic. 

Their result showed 
- 1/4 

p b t = constant 

This is actually a special case of Equation (18), it must 
be valid for any form of background field. 

Parker (1965), Belcher (1971), and Hollweg (1973b) studied 
the propagation of Alfven waves in the solar-wind with 
the simplification that the sun does not rotate so that 
the background magnetic field is in radial direction 
everywhere. They also assumed the medium to be. thermally 
isotropic and limited the fluctuation vectors to a linear 
polarization in the ^“direction of a spherical polar- 
-coordinate system. They obtained 

p '*' /,1+ (^±l)b^ = constant. 

Under the same assumptions, Hollweg (1973a) has also studied 
the higher-order terms, b 2 ,b 2 ,U 2 and u^, of the WKB approxi- 
mation. 

The accuracy of the WKB solutions can be gauged by comparing 
the magnitudes of the successive terms in the series repre- 
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senting the amplitude vectors. If we substract (15) from 
(14), we obtain 

“2-1*- 7 f 1 w ^2-l' 7f + fA-l - ^^ + M*V±Ae x . v)gp 


= 0 


(24) 




If and b 2 are parallel vectors, they must be parallel 
to u-^ and b^. The relative magnitude of the three terms 
in (24) may be measured by 




u- 



. 

■ 8n A c° s <f, 


+ 



The amplitude vectors may be approximated by their first 
terms , b^ and u^ , in the domain where 

pb 2 /bi < l 

and 

WUj/u^l 

These conditions are satisfied only when 


i-(_U 

8tt A cos$ 


+ 


Dl<l 


(25) 


The relation between the two fluctuation vectors described 
by Equation (8) holds under this approximation. In the 
region outside 1 a.u. the left hand side of (25) remains 
almost independent of the heliocentric distance r, because 
in that region A cos<(> is proportional to r ^ and 

In general the microscale fluctuation vectors may be re- 
presented by series of monochromatic waves. The intensity 
of the microscale fluctuations may be defined as the root- 
-mean-square value 


b 


( <b 2 > ) 


1/2 


(26) 
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Now we use the WKB solution. Equation (23), to calculate 
the variation of b in the interplanetary space. A self- 
-consistent solution of the solar-wind model including 
Alfvcn .waves and spiral fields is not available at the 
present time. In order to estimate the variations of £, p , 

U, A, and <p as functions of heliocentric distances, we 
choose the quantitative solution of the two-region model 
of the solar wind (Whang, 1972) in which the plasma thermal 
anisotropy and the spiral background field have been in- 
cluded. Figure 2 shows that b is calculated as a function 
of r based on the assumption that 

b = 3 gammas at 1 a.u. 

b varies from 370 gammas at 0.01 a.u., 50 gammas at 0.1 a.u. 
to 0.1 gammas at 10 a.u. 

The relative intensity of b in comparison with the magnitude 
of the background field, -B , is a more important parameter 
in describing the local condition of the solar wind. From 
Figure 2 we can see that the relative maximum of b/B Q is 
located in the region near 1 a.u. Thus outside of this 

region f for the Alfven mode in concerned, the solar wind is 

expected to be less turbulent. The relative intensity b/B 

- 1/2 . . ° 
varies asymptotically as r at large hiliocentric 

distances. 

IV GOVERNING EQUATIONS 

The equation of motion which governs dynamics of an 
anisotropic medium in the micro-, meso- and macro-scale is 

a 5 

~(pV) + V* (pVV + P I + — I - — — B B) - pg = 0 (27) 

at - k — = g, tt = 4* — 

Making use of (1) and (2), we can write (27) as 

!^(pU)+v- (p V U + pu u + 

+ p I + THT I - fcl ®o - & *!> - P§] * 0 <2S) 




7u - 


4ttP 


B • 7 (£b) 
— o — 


11 


No terms have been dropped out in Equation (28). The 
terms inside the first pair of square brackets is identical 
with Equation (3) which governs dynamics of microscale 
Alfven~ waves . Now taking the time average of Equation- (28), 
we obtain 

3U • B I P 

p( __ + U.VU) + 7.(P ± ! + l- 4^00^0+1 )-p£ = 0 (29) 

“r* 1 2 r 

where 1 1 = -5— < b >1 + p<u u> - 7*— <b b> 

= 8 71 = K 4 tt 


All terms in this equation are actually contributed from 
the time average of the terms inside the second pair of 
square brackets in (28). It governs the meso- and macro- 
scale motion of the background flow. In writing (3) as a 
governing equation for the microscale motion superposed 
on the background flow, no assumptions regarding the small- 
ness of the fluctuation vectors have been introduced. The 
stress tensor T is the only coupling term between the micro 
scale fluctuations and the background flow. To the zeroth- 
order approximation, this fluctuation stress corresponds 
to an additional isotropic^ressure due to microscale 
fluctuations (Alazraki and Couturier, 1971; Belcher, 1971) 


P A ' 8¥ <t>2> • 

Faraday's law in the micro-, 
be written as 


meso- and macroscale may 


3B 

TF + ' I!) = 0 


or 


[ft + '•W S - B^)] ♦ [■ 
- b u)] = 0 


8 B 

z2 +7.(V B - B U + u b - 

3 1 — * — O — — — — 


(30) 


The terms inside the first pair of square brackets is again 
identical with Equation (4) for the microscale Alfven waves 
The time average of the second part gives 
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3 ®o 

+ v - BqU ) = V • ( <b u> - <u b> ) 


(31) 






This equation governs the temporal and the spatial 
variations, of the background field in the meso- and. macro- 
-scale. To the zeroth-order approximation, the right hand 
side of (31) vanishes when the two fluctuation vectors, 
b and u, are either parallel or anti-parallel. 

It is now clear that when we write (3) and (4) to describe 
the microscale Alfven waves superposed on a background 
flow governed by (29) and (31), no limitations on the 
magnitude of the fluctuation vectors have been imposed. 

The WKB solution obtain in this paper, which satisfies 
(3) and (4) in a approximate way, must be valid for Alfven 
waves with no limitation on the amplitude. This special 
feature of the Alfven waves (Alfven and Falthammar, 1963) 
is preserved throughout this paper. 

To the zeroth-order approximation, the time-average 
energy equation, which governs the background flow in 
meso- or macro-scale, is 

' 2 r. 

+ £) + v. < 

e 




(|u 2 + 2P ± + 2* l + 


B 

4¥>U - 


-B B • U + q + S 
4tt— o— o — - 1 — 


- p j* • U - 0 


(32) 


It contains two additional terms contributed by microscale 
fluctuations. They are the fluctuation energy density 


£ = W> + 


8tt 


1 <b 2 > 


(33) 


and the total energy flux due to fluctuations 


S =SU + P^U - Ap <u > 


(34) 


The total flux consists of three parts: convection of the 
fluctuation energy density, rate of work done by the fluc- 
tuation pressure, and propagation of wave energy along the 
moving background field lines. Recently Hollweg (1973a) 



theorized that very low-frequency Alfven waves carry a 
larger wave energy flux than would be obtained in the 
zeroth-order WKB limit. 

In> the interplanetary* space, the energy flow due to- micro- 
scale Alfven waves per steradian per unit time, r S , is 
a monotonically decreasing function of hiliocentric 
distance. Its value near the corona ( 0.01 a.u.) is 
about 14% of the convective kinetic energy flow per 
steradian at 1 a.u. when we assume that the fluctuation 
intensity b = 3 gammas at 1 a.u. The estimated percentage 
of energy flow rises to 24% for b = 4 gammas at 1 a.u. 

This energy is converted into kinetic energy to increase 
the solar wind velocity. The interaction between the solar 
-wind flow and the attenuation of the fluctuation inten- 
sity is nonlinear in nature. A self-consistent model of 
the solar-wind must be studied to evaluate the effect 
of their interaction on -a more accurate quantitative basis 
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FIGURE. CAPTIONS 

Figure 1 Propagation of Alfven waves in a convecting 
spiral interplanetary magnetic field. 

Figure 2 The intensity of field fluctuation due to 
microscale Alfven waves is calculated as a 
function of heliocentric distance based on the 
assumed 3 gammas fluctuation at 1 a.u. The 
dashed curve is the magnitude of the background 
magnetic field. 
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